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ABSTRACT 

We analyze in some detail a recently proposed transfer matrix mean field ap- 
proximation which yields the exact critical point for several two dimensional nearest 
neighbor Ising models. For the square lattice model we show explicitly that this 
approximation yields not only the exact critical point, but also the exact boundary 
magnetization of a semi-infinite Ising model, independent of the size of the strips 
used. Then we develop a new mean field renormalization group strategy based on 
this approximation and make connections with finite size scaling. Applying our 
strategy to the quadratic Ising and three-state Potts models we obtain results for 
the critical exponents which are in excellent agreement with the exact ones. In this 
way we also clarify some advantages and limitations of the mean field renormaliza- 
tion group approach. 



1. Introduction 

A recently proposed [1] transfer matrix version of a mean field approxima- 
tion (which in the following will be denoted by LS) applied to several nearest 
neighbor Ising models in two dimensions, gave surprisingly exact results for 
the critical points, even without extrapolation, and very good results, under 
extrapolation, for more complicated models. 

A first issue to address in connection with the LS approximation is why 
the results are exact in the n.n. Ising case and extrapolate accurately in the 
others, and whether extra exact results can be obtained by this scheme. 

We want also to clarify what is the connection of this method with other 
techniques of more common use in two dimensional statistical mechanics. In 
particular: 

a) since the method involves consideration of strips similar to those used, 
e.g., in finite size scaling (FSS) and phenomenological renormalization 
approaches [2], it is legitimate to ask up to which extent the LS approxi- 
mation is connected to these approaches and possibly fits within them; 

b) since the method uses as a basic ingredient effective fields on the boundary 
of the strips, it is also rather natural to look for connections with the so- 
called mean field renormalization group [3] (MFRG) approach. This will 
go together with showing how critical exponents can be obtained in this 
context. 

The plan of the paper is as follows: in Sec. 2 we give a brief review of 
the LS approximation and show, in the square lattice Ising case, that it gives 
not only the exact critical point, but also the exact boundary magnetization 
of a semi-infinite Ising model, independent of the size of the strips used; 
in Sec. 3 we show how the LS approximation fits into a MFRG structure, 
develop a procedure for calculating the critical exponents, and compare the 
method with the FSS approach; in Sec. 4 we give two test applications of the 
procedure above, on the two dimensional Ising and three-state Potts models. 
Finally, in Sec. 5, we draw some conclusions. 
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2. The LS Approximation 

The LS approximation scheme [1] makes use of two infinite strips Sn and 
Sn' of widths n and n' respectively, with periodic bomidary conditions along 
the infinite direction. The approximation is obtained by applying an effective 
field hef f at one side of the strips and by imposing the consistency relation 



where rriin and niin' are the values of the order parameter at the opposite 
side of the strips and K = (3 J > is the exchange interaction strength. Eq. 
(1) has to be solved for he/f at fixed K and the critical temperature is the 
one for which the paramagnetic solution hgff — bifurcates into nonzero 
solutions, leading to spontaneous magnetization. 

This method, as pointed out by Lipowski and Suzuki [1], yields the exact 
critical temperature of the Ising model with nearest-neighbor interaction on 
many two dimensional lattices (square, triangular, honeycomb and centered 
square), and very accurate estimates of the critical temperature of more com- 
plicated models (Ising model with alternating strength of interaction, with 
next-nearest-neighbor interaction, 5' > 1 models). 

In the present section we will show, resorting to a result by Au-Yang 
and Fisher[4] , that at least in the simplest case of the nearest-neighbor Ising 
model on the square lattice the LS approximation yields not only the exact 
critical temperature, but also the exact boundary magnetization of the semi- 
infinite model. 

Let us consider an Ising model on a strip Sn of a square lattice, number 
from 1 to n the chains which form the strip and apply a magnetic field h^, on 
the nth chain. The corresponding hamiltonian will be 



where Sij = ±1 is an Ising spin located at the site with coordinates i and j in 
the X and y direction respectively. The magnetization niin = min{K, hn) = 



min{K, heff) = min'{K, heff), 



(1) 




(2) 



2 



(sj^i) has been calculated for n > 2 by Fisher and Au-Yang in ref. [4], and 
is given by 



min{K,hn) = z' {^-^ 



1/2 



|t| - t tanh(2n sinh"^|t'|) ^^'"^ 



■ (3) 



\i\ + (t + c+z'^)tanh(2n sinh ^\t'\) 
In eq. (3) we have adopted the same notation of ref. [4] , i.e. 

z' = taxihhn, 

= (^^i^' f= +«'^)'/^ = im -cosh2i^). 

c± = (2 + + t'^) -i±l^ cosh2K ± 1. 

For K > Kc = + -\/2) it is easily realised that if hn is chosen in such a 

way that t + c+z' = —t, i.e. 

, /cosh2K -coth2K\^/^ 

z' = tanhhn = ; (4) 

V cosh2K + l J ^ ' 

then the quantity in square brackets in eq. (3) equals 1, independent of n, 
and one has 

/ cosh2K-coth2i^ V/' 
mUK.K) = m,{K) = eosh2K - 1 ) 

where mi{K) is the exact boundary magnetization of the two dimensional 
semi-infinite Ising model[5]. Furthermore, for K < Kc, choosing hn = 
yields mi{K) = 0, again independent on n. 

The case n = 1, for which eq. (3) is no more valid, can be easily solved 
by the transfer matrix method, obtaining 

mii{K, hi) = (6) 

(sinh^/ti + ^ 

which again equals mi{K) if hi is chosen according to eq. (4). 

These results imply that, no matter which n, n' we choose, the bifurcation 
at heff = will always occur at the exact critical point, K = Kc- This 
mechanism is at the basis of the results of ref. [1] for the square lattice case. 
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and we believe that the same should work for the other two dimensional 
lattices. 

Finally, if, fully in the spirit of a classical approach, we consider the 
nonzero magnetization solution in eq. (1) for K > Kc, by putting heff = hn 
we obtain the exact spontaneous boundary magnetization mi, independent 
of n, n'. This was not noticed in ref. [1]. 

The following remarks are in order for an explanation of the above results. 
First of all, the boundary magnetization is known to behave as mi{K) ^ 
(K — Kc)^^ for K with (3i = 1/2, in the 2d Ising model. The exponent 

Pi = 1/2 [5] is such that it can be reproduced exactly by a mean field like 
approximation, making use of an effective field. When considering models 
with I3i values incompatible with a classical scheme, one has to consider 
the LS approach and its possible extensions and approximations, as we will 
discuss in the next sections. 

As shown in ref. [6], in the context of a generalized cluster variation 
approach to two dimensional lattice models, the double strip ^2 is able to 
contain the whole information needed to solve exactly the two dimensional 
n.n. Ising model. The problem then reduces to how such information can be 
extracted. Clearly what we presented in this section amounts to a relatively 
simple way of obtaining part of this information. 

3. Mean field renormalization group and finite size scaling 

Let us now see how the LS approximation can be used to develop a new 
MFRG strategy, where the boundary magnetization is used together with 
the bulk one as an effective scaling operator. This will also be useful in 
understanding the relations of the LS method with FSS. 

The notation applies to an Ising model for convenience, but the strategy 
is not limited to this will be shown in the next section where it will 

be used to investigate the three-state Potts model. 

For infinite Ising strips of widths n and n', FSS implies the following 
scaling law for the singular part of the bulk free energy density, f^^\ 

fi^\£y-e,£y-h)^£''fi'\e,h), (7) 
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where e = (Tg — T)/Tc., i = n/n' is the rescahng factor, and d is the bulk 
dimension. If the boundary conditions are open, for the singular part of the 
surface free energy density, f^^\ the relation 

fi',\ey-ejy"h,ey"'b) = fi'\e,h,b), (8) 

holds, where b indicates a surface field. 

The MFRG basic idea[3] is to derive from eq. (7) the scaling relation for 
the bulk magnetization 

mn'{K',h') = e''-y"mn{K,h), (9) 

where h' = i^" h and K' = K'{K) is a mapping in the Wilson-Kadanoff sense, 

determined implicitly, in the limit of h going to zero, on the basis of eq. (9). 

From this mapping the critical point Kc and the thermal exponent j/y can be 

OK' 

obtained by means of the relations K(. = K'{Kc) and = = Kc). 

oK 

Applying this idea to the surface magnetization yields 

min'{K',h',b')=e''-^-y"'min{K,h,b), (10) 

where b' = i^^^b if we want b to scale as a surface field. On the other 
hand, the equation for the critical point which is characteristic of the LS 
approximation would be recovered if b scaled as a magnetization, i.e. with 
an exponent d — 1 — yns- In fact, with this assumption, setting /j, = and 
linearizing in b yields 

^(jr',o,o) = ^(if,o.o), (11) 

which implicitly defines a mapping K\K). The equation 

Kc = K\Kc) (12) 

with K'{K) given by eq. (11) is equivalent to the equation for the critical 
point in the LS approximation. So this approximation can also be seen as a 
realization of a MFRG strategy as far as determination of Kc is concerned. 
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In a MFRG spirit one can also determine the critical exponents, since j/y 
is obtained by the relation 

dK' 



dK 



(13) 



linearizing eq. (9) with respect to /i, with K = Kc yields 

^(K„0)=£'^-2--^(Ke,0), (14) 

from which hh can be obtained and finally, linearizing eq. (10) in the same 
way, with 6 = 0, yields 

^(Ke, 0, 0) = ^'^-^-VHs-VH^J^^K,^ 0, 0), (15) 

from which ijhs is obtained. 

The set of equations (11)-(15) is a MFRG procedure to determine critical 
point and critical exponents. 

Nevertheless, the procedure above (to be denoted by M, for "magneti- 
zation", in the following) is not a rigorous application of FSS. In such an 
application (11) should be replaced by 

^(K', 0, 0) = ^(^, 0, 0), (16) 

since h should scale as a field, with exponent yns-, and should be solved in 
conjunction with (14)-(15). This alternative and more rigorous procedure 
will be denoted by F, for "field", in the following. 

F is in fact the procedure of MFRG proposed in ref. [3] to yield simulta- 
neously bulk and surface exponents. It is interesting to investigate how M, 
proposed here, being more consistent with the effective field idea, compares 
with F. 

Two comments are in order: 
i) the two procedures should give the same value of K^. (but not of the 

exponents) in the limit n, n' ^ oo, £ ^ 1, since the two derivatives in eq. 

(11) are analytic functions; this should justify the LS approximation in a 

FSS context; 
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ii) in the two dimensional n.n. Ising case, yns = 1/2 exactly and then, in 
the limit above, also the critical exponents should be the same for both 
procedures M and F. 

In the following section we will check these ideas on the two dimensional 
n.n. Ising and three-state Potts cases. 



4. Results and discussion 

In the present section we give two test applications of our new MFRG 
strategy (M), to the Ising and three-state Potts models on square lattices. 
We also compare our results with those obtained treating 6 as a surface field, 
i.e. letting it scale with exponent yns- We start by applying procedure M to 
the Ising model. In the Ising case, we have already shown that the method 
gives the exact critical point for any n, n' . Furthermore, resorting to eq. (3), 
the mapping K' = K'{K) can be determined analytically in an implicit form. 
As a result one gets 

fn'{K') = U{K), (17) 



where 

1/2 

I -J- T^nnt /Ti Ginn it i i 

fn{K) = K coihK 



1 + tanh(2n sinh ^\t'\) 
1 — tanh(2n sinh""^ l^'l). 



(18) 



and with t' as above. It can be checked that the fixed point of (17), obtained 
by setting t' = 0, is K* = ^ ln(l + a/2), while for the thermal exponent one 
has 

^ l + 2(2n-l)K* 

l + 2(2n'-l)K*' ^ ^ 

which, in the limit n,n' — > oo, yields j/t = 1, which is again an exact result. 

The calculation of the magnetic exponents cannot be carried out analyt- 
ically since no solution is available for the bulk magnetization of a strip in 
the presence of a bulk magnetic field, and we have to proceed numerically, 
as follows. Given the strip Sn with a bulk magnetic field h and an auxiliary 
magnetic field hi acting on the first chain, we determine its partition func- 
tion Zn{K, h, hi) as the largest eigenvalue of the 2"^ x 2"^ transfer matrix with 
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elements 



Tn{{sj}As'j}) = e^p 



K 



n-l 



h 



hi 



(20) 



The bulk and boundary magnetizations will then be given by 

1 dZn 



and 



rUn = mn{K, h) 



min = min{K, h) 



Zr,. dh 



1 dZr, 



(21) 



(22) 



/n=o 



Z„ dhi 

respectively. Finally yn and yns are determined according to (14)-(15). 

In Tab. 1 we report the results for the critical exponents for strip widths 
2<n = n' + l<8(n = n' + lis always the most convenient choice). The 
extrapolations are based on least squares fits with fourth order polynomials in 
1/n and are certainly justified since the critical exponents are nearly linear 
functions of 1/n. The agreement of the extrapolated data with the exact 
results is very good, and the errors are within 0.2%. 

The Ising test has yielded very promising results, but does not shed much 
light on the physical meaning of the effective parameter 6, since for the two 
dimensional Ising model surface magnetization and surface field scale with 
the same exponent d — 1 — yns = Uhs = 1/2. In fact, procedure F, in which 
b scales as a surface field yields as well very good results (apart from having 
no solution when n = 2), as shown in Tab. 2. Data from the two procedures 
are plotted together in Figs. 1-4. As expected, all the results seem to be 
equivalent in the limit n — > oo. 

In view of the above considerations, we believe that a more conclusive test 
is in order, and a suitable model should be the three-state Potts model. In- 
deed, in two dimensions, this model is known to undergo a second order phase 
transition, whose critical point and critical exponents are known exactly [7,8], 
even in the absence of a full solution. In particular it has yns 7^ 1/2. 
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The hamiltonian of the q-state Potts model[7] is 



K 



1) + 



h 



(23) 



q-1 



q-1 



where the variables Si take on values 0, — 1, K > is the interac- 

tion strength and h is a magnetic field. The order parameter of the model, 
corresponding to the Ising magnetization, is 



In the case q = 2 one recovers the Ising model. 

The MFRG scheme developed above can be carried over to the ^--state 
Potts model without any substantial modification, and we will apply it to 
the case q = 3. The main new fact is that no analytical results like eq. (3) 
are available for the three-state Potts model. So all calculations must be 
performed numerically with the transfer matrix method. The order of the 
transfer matrix is now 3"^ and increases more rapidly than in the Ising case. 
However, the transfer matrix is invariant with respect to the transformation 
which interchanges the states Si = 1 and Si = 2 (all other symmetries are 
lost as soon as one introduces the surface fields), and the eigenvector corre- 
sponding to its largest eigenvalue belongs to the symmetric subspace of this 
transformation. Thus we can limit ourselves to matrices acting in this sTib- 
space, which are of order (3"^ -|- l)/2. In this way we have been able to deal 
with strips up to n = 7. 

The numerical results for procedure M are reported in Tab. 3. The ex- 
trapolation is obtained by fitting data in a least square sense with a second 
order polynomial in 1/n. Even if now the critical point is not given exactly 
by the LS approximation, we obtain excellent agreement with the exact re- 
sults (error within 2.3%) for the bulk critical point and exponents, while, in 
comparison, the results for the surface exponent i/hs are rather poor. 

In Tab. 4 we report the results obtained from procedure F. The situation 
is very similar to the previous one: the errors of the extrapolated critical 
point and bulk exponents are within 4.7% and, in comparison, the estimate 



m = 



q-1 



(24) 
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of the surface exponent is again rather poor. Results from the two procedures 
are plotted together in Figs. 5-8. 

There is some evidence that F is the correct procedure when n is large, 
but for small strips M, although not rigorous, seems to work very well: indeed 
if i/H had been extrapolated on the basis of the results for 2 < n < 5 one 
would have obtained 1.868, which is two orders of magnitude more accurate 
than the extrapolation on the whole set of data. 

5. Conclusions 

We have analyzed in some detail the LS approximation, showing that in 
the two dimensional n.n. Ising case, it yields not only the exact critical point, 
but also the exact boundary magnetization of the semi-infinite model, inde- 
pendent of the size of the strips used. We have also proposed an explanation 
of these surprising results. 

The LS approximation has been used to develop a new MFRG strategy 
(procedure M) which yielded very accurate results for the critical exponents 
of the Ising and three-state Potts models in two dimensions. When com- 
pared with rigorous FSS (procedure F) our new strategy has proven to be 
particularly suitable for applications where relatively small strips are used. 
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Figure Captions 



Fig. 1 : The Ising critical point Kc vs. 1/n as given by the M (sohd hne) 

and F (dashed hne) procedures. 
Fig. 2 : The Ising thermal exponent ut vs. 1/n as given by the M (solid 

line) and F (dashed line) procedures. 
Fig. 3 : The Ising magnetic exponent yn vs. 1/n as given by the M (solid 

line) and F (dashed line) procedures. 
Fig. 4 : The Ising surface magnetic exponent yns vs. 1 /n as given by the 

M (solid line) and F (dashed line) procedures. 
Fig. 5 : The three-state Potts critical point K(. vs. 1 /n as given by the M 

(solid line) and F (dashed line) procedures. 
Fig. 6 : The three-state Potts thermal exponent j/t vs. 1/n as given by 

the M (solid line) and F (dashed line) procedures. 
Fig. 7 : The three-state Potts magnetic exponent yu vs. 1/n as given by 

the M (solid line) and F (dashed line) procedures. 
Fig. 8 : The three-state Potts surface magnetic exponent yns vs. 1/n as 

given by the M (solid line) and F (dashed line) procedures. 
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Tab. 1: results for the Ising model (6= magnetization) 

n Vt Vh Vhs 



2 


0, 


.95738 


1 


.60287 


0. 


,60287 


3 


0, 


.97310 


1 


.68092 


0, 


,56808 


4 


0, 


. 98089 


1, 


.72197 


0, 


.55099 


c 
O 


U . 


. yt5oi4 


1 , 




U . 


. t)4UD ( 


6 


0, 


. 98785 


1, 


.76600 


0, 


.53374 


7 


0, 


.98969 


1, 


.77937 


0, 


.52873 


8 


0, 


.99106 


1, 


.78970 


0, 


, 52496 


Extrap 


0, 


.99977 


1, 


.87201 


0, 


, 49764 


Exact 


1 




1, 


.875 


0, 


,5 



Tab. 2: results for the Ising model (b — field) 



n 




Kc 




VT 




VH 




Vhs 


3 


0, 


.52903 


0, 


. 96548 


1, 


.99012 


0, 


.91150 


4 


0, 


.48065 


0, 


. 97840 


1, 


.91195 


0, 


.77040 


5 


0, 


.46366 


0, 


. 98404 


1, 


.88570 


0, 


.70267 


6 


0, 


.45560 


0, 


. 98726 


1, 


.87397 


0, 


.66195 


7 


0, 


.45113 


0, 


. 98938 


1, 


. 86794 


0, 


. 63455 


8 


0, 


. 44839 


0, 


.99086 


1, 


.86461 


0, 


.61482 


Extrap 


0. 


.44512 


0. 


.99858 


1, 


. 87440 


0. 


.49865 



Exact . 44069 1 



1 . 875 0.5 



Tab. 3: results for the 3 — state Potts model (b = 
magnetization) 



n 




Kc 




Vt 




Vh 




Vhs 


2 


0. 


,71318 


1, 


.10152 


1, 


.60916 


0, 


,60916 


3 


0, 


.69762 


1. 


,12118 


1, 


.69079 


0, 


.56966 


4 


0, 


.69013 


1. 


,13209 


1, 


.73652 


0, 


. 54785 


o 


n 

u . 


. DOO f '± 


1 . 


, ±oy I o 


i 






. OooOU 


6 


0, 


. 68287 


1. 


, 14496 


1, 


.78923 


0, 


.52317 


7 


0, 


. 68084 


1. 


, 14993 


1, 


.80629 


0, 


.51526 


Extrap 


0, 


.66896 


1. 


,17561 


1, 


.90880 


0, 


.47250 


Exact 


0, 


. 67004 


1. 


,2 


1, 


.86667 


0, 


.33333 



Tab. 4: results for the 3 — state Potts model (b — field) 



n 




Kc 




VT 




VH 




Vhs 


3 


0, 


.78999 


1. 


,12871 


1, 


.96588 


0, 


.87340 


4 


0, 


.72612 


1. 


,14060 


1, 


.88767 


0, 


.72035 


5 


0, 


.70302 


1. 


, 14627 


1, 


.86166 


0, 


.64290 


6 


0, 


.69184 


1. 


,15025 


1, 


.85030 


0, 


.59417 


7 


0, 


. 68555 


1. 


,15541 


1, 


. 84477 


0, 


.56018 


Extrap 


0, 


.70167 


1. 


,17551 


1, 


.88526 


0. 


.40703 



Exact 0.67004 1.2 1.86667 0.33333 



